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Abstract 

In this work, we explore some cosmological implications of the model 
proposed by M. Visser in 1998. In his approach, Visser intends to take 
in account mass for the graviton by means of an additional bimetric tensor 
in the Einstein's field equations. Our study has shown that a consistent 
cosmological model arises from Visser's approach. The most interesting 
feature is that an accelerated expansion phase naturally emerges from the 
cosmological model, and we do not need to postulate any kind of dark energy 
to explain the current observational data for distant type la supernovae 
(SNIa). 
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1. Introduction 

The state of the art in cosmology has led to the following present 
distribution of the energy densities of the Universe: 4% for baryonic matter, 
23 % for non-baryonic dark matter and 73 % for the so-called dark energy 
(see e.g. pQ). The dark components (matter and energy) have been the 
focus of studies and a lot of speculations in theoretical and observational 
astrophysics. 

In particular, the dark energy is very curious, not only for its dominant 
relative energy density, but for its dynamical consequences on cosmic 
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expansion: it can make the Universe accelerate. Phenomenologically, some 
equations of state have been proposed in order to explain such a dark energy. 
The most common example is the cosmological constant (ACDM model), 
which implies on a constant vacuum energy density along the history of the 
Universe. However, the quantum theory prediction for the vacuum density 
is 120 orders of magnitude greater than the observational value (see e.g. [2]). 

Moreover, ACDM model is not satisfactory because it requires a large 
amount of fine tuning to produce cosmological constant energy density 
dominant at recent epochs (see e.g. [3]). 

Another possibility is a dynamical vacuum or quintessence. In general, 
the quintessence models involve one [I] or two [5] coupled scalar fields. 
In these models, the cosmic coincidence, i.e., why dark energy started 
dominating the cosmic evolution only so recently, has no satisfactory solution 
and some fine tuning is required. 

The Chaplygin gas is another example of dark energy fluid. Its exotic 
equation of state can be derived from the Nambu-Goto action for 'D-Branes' 
moving on a (d+2) dimensional space-time [6]. The Chaplygin Gas has a 
dual behavior: in the past it behaves like matter and in recent times like 
a cosmological constant. So, it can represent dark matter and dark energy 
with only one equation of state. However, at the time of structure formation 
the influence of dark energy component is negligible and matter clustering 
occurs in a similar way as in CDM model [7]. 

More recently some authors (see, e.g., [3]) have studied the dark energy 
problem considering a quintessence model with dark matter - dark energy 
interaction. However, as discussed by Bertolami et al [8], such a kind of 
coupling between dark matter - dark energy produces a violation of the 
equivalence principle (EP). Actually, a violation of the EP is reported to 
be found in other dark energy models [9]. 

Whatever the dark energy may be, it seems that physics beyond the 
standard model is necessary. 

On the other hand, several studies have been trying to bound the mass 
of the graviton by many ways. From analysis of the planetary motions in 
the solar system it was found that we must have m g < 7.8 x 10~ 55 g pU] 
in order to respect the accuracy of the observations with the newtonian 
potential. Another bound comes from the studies of galaxy clusters, which 
gives m g < 2 x 10~ 62 g [llj. Although this second limit is more restrictive, it 
is considered less robust due to uncertainties in the content of the Universe 
in large scales. 
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Studying rotation curves of galactic disks, de Araujo and Miranda [12] 
have found that m g <C 10~ 59 g in order to obtain a galactic disk with a scale 
length of b ~ 10 kpc. 

Studying the mass of the graviton in the weak field regime Finn and 
Sutton have shown that the emission of gravitational radiation does not 
exclude a possible non null rest mass. They found the limit m g < 1.4 x 10~ 52 g 
[T3] analyzing the data from the orbital decay of the binary pulsars PSR 
B1913+16 (Hulse-Taylor pulsar) and PSR B1534+12. 

As can be seen, the graviton mass is not observationally excluded. So, it 
is reasonable to ask if the consideration of mass for the graviton plays some 
role in cosmology. Do cosmological models exclude a non null graviton mass? 
Can massive gravitons affect the cosmic dynamics? Is it possible to bound 
the graviton mass by cosmological observation? Would such a bound be in 
accordance with other observations? Considering mass for the graviton, do 
we need to include dark energy in order to explain the cosmological data? 

In order to look for answers to these questions, we have chosen to include 
the graviton mass adopting the Visser's approach [H], where the graviton 
mass appears as an extra term in the Einstein's equations. It is worth noting 
that the weak field equations used by Finn and Sutton are the same that 
come from the Visser's model when we use the linear approach. 

An interesting result that comes out from this model is that the 
gravitational waves present six polarization modes [TJ)] instead of the two 
usual polarizations obtained from the General Relativity theory. So, if in the 
future we would be able to identify the gravitational wave polarizations, we 
would impose limits on the graviton mass by this way. 

This paper is organized as follows: in Section 2 we briefly review the 
Visser's approach. Section 3 is devoted to the description of the cosmological 
model. In Section 4 we discuss how the age of the Universe can constrain 
the value of the mass for the graviton. Section 5 describes the evolution 
of the scale factor in the 'massive cosmology'. In Section 6 we show that 
a cosmological scenario with massive gravitons can produce a phase of 
accelerating expansion for the Universe. Section 7 presents a comparison 
between ACDM and our cosmological model using the luminosity distance. 
Finally, in Section 8 we present our conclusions. 
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2. The Field Equations 

The full action considered by Visser is given by [Hj: 

-c A R(g) 



d x 



-g- 
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where besides the Einstein-Hilbert lagrangian and the lagrangian of the 
matter fields we have 
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where m g and {go)^u are respectively the graviton mass and a general flat 
metric. 

The field equations, which are obtained by variation of can be written 

as: 



CT - \ = -^T^, (3) 

where G^ v is the Einstein tensor, T^ v is the energy-momentum tensor for 
perfect fluid, and the contribution of the massive tensor to the field equations 
reads: 



= (g^r [(g - g U - ^oMsq- 1 )^ 
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Note that if one takes the limit m g — > the usual Einstein field equations 
are recovered. 

Regarding the energy-momentum conservation we will follow the same 
approach of [TU] and [T7] in such a way that the conservation equation now 
reads [181: 



i f> 
m g c 



since the Einstein tensor satisfies the Bianchi identities V V G^ U = 0. 



(5) 
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3. Cosmology with massive gravitons and without dark energy 



For convention we use the Robertson- Walker metric as the dynamical 
metric: 

^ r i 2/i/i2 i • 2qj±2\ 

— + r (at) + sin " 1 1 



ds 2 = c 2 dt 2 - a 2 (t) 



(6) 



where a(t) is the scale factor. The flat metric is written in spherical polar 
coordinates: 

ds 2 = c 2 dt 2 - [dr 2 + r 2 (d6 2 + sin 2 Odcf 1 )] . (7) 

Using ^ and ^ in the field equations ^ we get the following equations 
describing the dynamics of the scale factor (taking k = for simplicity): 
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where as usual p is the energy density and p is the pressure. 

From equation ^ we get the evolution equation for the cosmological 
fluid, namely: 



P + 3H 
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where H = a/a. Considering a matter dominated universe (p = 0) the above 
equation gives the following evolution for the energy density: 
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where p an d a are the present values of the energy density and the scale 
factor respectively. Note that in the case m g 4 we obtain the usual 
Friedmann equations. 

Now, inserting (11) in the modified Friedmann equation ^ we obtain 
the Hubble parameter: 



H 2 {a) = Hi 



(12) 
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where the relative energy density of the 2-component is Q{ = Pi/p c (p c = 
3H 2 c 2 /8irG is the critical density) where 'i' applies for baryonic and dark 
matter. Here ran = hHo/c 2 is a constant with units of mass, which we will 
call 'Hubble mass' (as we will see later this constant is very important in the 
present context). 



4. The Age of the Universe Constraining the Mass for the Graviton 



The equation (12) for the present time (a = a ) gives: 

2 



£1° + 



( 7a g- 5a $) = l. 



140 \m H 

Solving this equation to a = a (m g ) we get: 



(13) 
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(14) 



Thus, if we use (14) in the equation (12) we have the dynamical equations 
described in terms of the free parameters f^, H and m g . 



Note that in order to have real values for ao, the term in brackets in (14) 
must satisfy the relation: 



which lead us to a lower limit for the graviton mass in our model: 



m g > 



200 



(1 - noj m H . 



If we take, for example, Vt° m = 0.27 we have: 

m g > 4.57mjy. 



(15) 



(16) 



(17) 



So, the Hubble mass, m H , establishes the magnitude of the expected value 
of the graviton mass. If we convert this lower limit in a upper limit to the 
Compton wavelength of the graviton we have: 



X g < 0.22 — 



(18) 
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Figure 1: Relation between the age of the Universe and the mass for the graviton for 
different values of Q, m . 

Thus, if the graviton rest mass is non null, its corresponding Compton 
wavelength must be lower than the observable horizon. But once the 
Compton wavelength is associated with the range of the interaction, this 
tell us that the contribution of the mass term can be relevant in the Universe 
at the present time. 

A note about the lower bound on the graviton mass. Since we are 
considering only matter in the context of a massive gravity theory, the lower 
bound on the graviton mass is a direct consequence of the fact that we have 
fltotai — 1 today. Once matter contributes only with 27% of the total relative 
energy density, and the radiation contribution is negligible, we could not 
have a zero graviton mass without considering another kind of dominating 
component, but such idea is not the aim of the present paper. 

Once we have an analytical description of the dynamics in this scenario 
we can impose limits on the parameters in many ways. One of them can be 
the age of the Universe. 

The time scale to the age of the Universe is related to the Hubble 
parameter, namely, tjj = Hq 1 , which is about 13 Gy as given by the current 
observational data pp. 

So, we hope that a consistent cosmological model can give ages of this 
order, which contemplates the age of the oldest stars and the time to structure 
formation. 
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In order to calculate the age of the Universe in our model we just solve 
the integral: 

tu = / -^-r- (19) 
J an [a) 



Substituting (12) in (19) we obtain the age of the Universe as shown in 
Figure [l]to some values of TO^. We can see that the age of the Universe is very 
closely related to the mass of the graviton. This gives us a very restrictive 
upper limit. The effect of taking different values of is to shift the curve 
upward or downward, but the qualitative behavior is the same in all cases. 

If we take, for example, tn = 13 Gy and = 0.27 we obtain the limit: 



m 



, g ho 1 < 2.24 x 1(T 65 g , (20) 

where the Hubble parameter is given by Hq = lOO/io kms~ l Mpc~ l and 
ho is a dimensionless constant, which parameterizes the uncertainty in the 
measurement of Hq. This limit for m g is about 10 orders of magnitude more 
restrictive than the limits obtained from the orbital motion in solar system 
and about 3 orders more restrictive than the inferences from galaxy clusters. 

Considering the lower limit obtained previously, the mass of the graviton 
in this model must be in the interval: 

1.73 x 10~ 65 g < m g h 1 < 2.24 x 10" 65 g , (21) 

for $7^ = 0.27. This upper limit depends on the value taken for the age of 
the Universe, while both upper and lower limits depend on the value of f2°. 



5. Past and Future 



One can determine how the Universe evolves by integrating equation ( 12 ), 
which yields a(t) as shown in Figure [2] 

In the past and in the present the curve a(t) behaves like the evolution of 
the Friedmann models with a cosmological constant. However, the future is 
drastically different. The massive term contributes in such a way that in the 
past it generates a repulsion, becoming attractive in the future and leading 
the Universe to the so called big crunch. Moreover, although the evolution 
is qualitatively the same for any value of the graviton mass, it is notable the 
strong dependence between the time of evolution and m g . This emphasizes 
again the restricted values to m g which we are considering. 
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In the Figure [3] we can see the evolution of the Hubble parameter. In 
the past the function behaves like the Friedmann models and in the future 
H(z) — > 0. This shows us that the massive contribution becomes important 
in the late time for the history of the Universe. In fact, if we take a — > 



in equation (12) we see that the massive term is negligible when compared 
to the radiation and matter contribution. That is very important in order 
to have no change in the expansion rate of the Universe, for example, in the 
nucleosynthesis era (a detailed study of this issue will appear elsewhere). 

When the scale factor reaches its maximum value (the turning point 
Ztum), the Hubble parameter goes to zero H(z turn ) = 0. The smaller m g 
is, the greater z iurn is. This indicates the relation between the turning point 
and the range of the gravity interaction, once the Compton wavelength is 
inversely proportional to the mass of the particle. 




Figure 2: Evolution of the normalized scale factor with time (in units of the Hubble time) 
for different values of m„ 



6. Accelerating Universe? 

In order to study the evolution of the second derivative of the scale factor 
we will use the definition of the dimensionless desaccelerating parameter: 

* s -*- (22) 
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Figure 3: Evolution of the Hubble parameter as a function of the redshift for different 



values of ', 



We take = 0.27. 



Applying this to the equations ^ and (|9| we get: 
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The evolution of the desaccelerating parameter as a function of the 
redshift is shown in Figure |4| 

Note that an accelerating expansion at the present time is closely related 
to the value of the graviton mass. In the past, as expected, the Universe has 
a desaccelerated phase. The fast growth of q(a) in the future is explained 
by the reversion in the expansion as discussed above. When the Universe 
is contracting, the values of q(a) are the same as in the expansion phase, 
which lead us to conclude that the Universe has two accelerating phases in 
this model: the first when it expands and the second when it contracts. 

In order to verify the dependence between the present value of the 
desaccelerating parameter (go) and the value of m g we just set a = 1 in 



the equation (23), giving: 



-n° + 

2 m 



m. 



3ag 

IT 



(24) 



The curve q x m g is shown in Figure |5] If we assume VP m = 0.27, for 
example, we have a present accelerating expansion form^ 1 < 1.80xl0" 65 g. 
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It is interesting to observe that this result is compatible with older 
Universes and greater lifetimes, which indicates that the most probable value 
for m g would be closer to the lower limit (see equation (21)). 
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Figure 4: Evolution of the desaccelerating parameter q(z) for different values of m g . We 
take n° n = 0.27 
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Figure 5: Relation between the current value of the desaccelerating parameter and the 
graviton mass for different values of f2 m . 
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7. The Luminosity Distance 



In order to test this cosmological model, we calculate the luminosity 
distance and compare it with the ACDM best fit of the type la supernovae. 
The luminosity distance to a flat Universe (k = 0) is given by (see, e.g. 

my- 

f z dz' 

d L (z) = (l + z)cH^J (25) 



where h(z) = H(z)/H and (1 + z) 



-l 



In the ACDM model the equation (25) gives us: 



iKCDM (1 + Z)C f Z dz' 



d MJDM = I = . (26) 

H Jo ^(1 + z'f + tt° A 



To the massive model we can combine (|12|) with (|25|) to obtain: 
d", 



jmass 



[l + z)cH- 1 [ rfz'vTlM, (27) 
Jo 



where: 



and 



n ^') - mW^m • (28) 



M(z') = (m g /m H ) 2 [4/U - (ag/10)(l + z'f] 



If we use the relation between ao and m g , given by (14), the luminosity 
distance in our model has the same number of free parameters as in the 
ACDM model. 

In Figure [6] we plot some values of the graviton mass for the limits given 



by (21 ). We also plot the curve for the ACDM model, which is the best fit for 
the most recent observational data from the recession of type la Supernovae, 
i.e., Q° m = 0.27 and = 0.73 [20J. The comparison between the model and 
the SNIa data can be seen in Figure [7], where we use the definition of the 
distance modulus: 

H = m - M = 5 log(di) - 5 + A , (29) 

where m and M are the apparent and absolute magnitude respectively and 
A is the absorption in magnitude. 
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Figure 7: Distance modulus as a function of the redshift to the massive model compared 
with the SNIa data from Astier et al. (2006). 

As can be seen, there are not significant differences in the curves di x z, 
when we change the graviton mass in that range we have considered. 
Furthermore, we can give a value to m g which describes exactly the same 
curve given by the ACDM model, this value is around nigh^ 1 = 1.95 x 10 -65 
g. From Figure ^ we see that for this value of the graviton mass we have a 
positive desaccelerating parameter in the present (z = 0), so we can fit the 
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observational data for SNIa without a present acceleration in the expansion 
of the Universe, although we have an inevitable accelerating expansion era 
in the near past. 

Again, we would like to emphasize that there are no significant changes in 
the curve d^x z when we use other values of m g , which are in the established 
limit. In such a way there are some possible values of m g , which give current 
acceleration and other that give desacceleration. For example, taking the 
values used in the Figure [7j the lower one (1.74 x 10~ 65 g) gives current 
desacceleration (q > 0), and the greater one (2.23 x 10~ 65 g) gives a current 
cosmic acceleration (q < 0). 

In view of the above results we may conclude that it is not possible to 
decide between the massive model and ACDM model only by the type la 
Supernovae observations, in particular for the small redshifts in which the 
Supernovae have been observed (z < 1). 

8. Conclusions 

We have shown that within the context of a classical gravity theory with 
massive gravitons, we can obtain a consistent cosmological model which has 
acceptable values for the age of the Universe, furthermore, it can fit the 
present cosmological SNIa data without any kind of dark energy. 

The possible values for the graviton mass in this theory are in accordance 
with all established limits. 

We identify that there is degeneracy between the ACDM model and the 
massive model when we analyze the luminosity distance versus redshift. Such 
a degeneracy could be removed, in principle, by other tests of the massive 
model such as structure formation or analyzing the power spectrum of the 
cosmic microwave background radiation (CMBR). 

In future works we will study in detail the nucleosynthesis era within the 
context of the massive model, as well as to provide a statistical approach to 
the analysis of the SNIa data within this model. 

Once the graviton mass is introduced via an effective tensor as given 
by (j3j), we hope that the primordial graviton production would be different 
from those models which consider general relativity. So, the future detection 
of primordial gravitational waves will provide a way to test this alternative 
gravity theory. 

It is worth stressing that Gabadadze and Gruzinov (2005) have analyzed 
the instabilities of the background and ghosts produced by massive gravitons 
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in 4D Minkowski space-time. They conclude that a natural way to account for 
a massive graviton on flat space is to invoke theories with extra dimensions. 
However, Visser's model is truly continuous with general relativity (GR) in 
the limit of vanishing graviton mass. Together with Pauli-Fierz (PF) massive 
term, Visser's theory is the simplest attempt to incorporate mass for the 
graviton in a ghost-free manner. Moreover, the van Dam-Veltman-Zakharov 
discontinuity (vDVZ) [22] present in the PF term can be circumvented in 
Visser's model by introducing a non-dynamical flat-background metric [23J . 

Thus, Visser's model can help us to understand the influence of the mass 
for the graviton in cosmology and maybe it could shed some light on the 
following question: what is the dark energy? 
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